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The concept of partnership of potentials is studied in detail and in particular the non- 
uniqueness due to the ambiguity in the election of the factorization energy and in the 
choice of the solution of certain Riccati equation. We generate new factorizations from 
old ones using invariance under parameter transformations. The theory is illustrated 
^ ' with some examples. 
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1. Introduction 

o : 

f^*) , The so-called Factorization Method plays an important role in the search for quan- 

tum systems for which the spectrum of the Hamiltonian operators is complct 
O |. known.EiBB It is closely related with the existence of an intertwining operator. 

with Supersyjnmetric Quantum MechanicsQ and Darboux transformations in this 
last context.0 Moreover, these techniques have important gene 
dimensional spacesJS to higher order factorization operators] 
class of systems with partial algebraization of the spectrum,E3 
Actually, most exactly solvable potentials can be obtained by making use of an 
appropriate intertwining operator transformation. 
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$_i ' Looking for a factorization of a given Hamiltonian amounts to find a constant 

d and a solution of a Riccati differential equation for the superpotential function. 
Once a solution has been found, a partner potential is defined such that it has 
almost the same spectrum as the original one. However, the point is that with 
a different choice for the solution of the mentioned Riccati equation, a different 
partner potential is obtained. This fact has been shown_to be very useful for the 
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lotentials, an idea due to MiclnikLj and later developed in 
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Moreover, the ambiguity in the definition of a partner potential of a given one 
is twofold: firstly, due to the choice of the factorization energy d, which is not 
unique in general, and then, it arises the ambiguity in the election of the solution of 
the corresponding Riccati equation. These two ambiguities are more or less known 
and are implicitly used or mentioned in many papers. However, we feci that it is 
worth having a new look at the subject in its own right because its understanding 
allows to interpret certain facts treated in the literature as consequences of this 
undetermination. In addition, these ideas have a great influence in the same concept 
of partnership when applied to the subclass of Shape Invariant partner potentials. 

Shape Invariance is an important concept in the theory of exactly solvable sys- 
tems, which was explicitly introduced by GendenshteintJ although the basic idea 
was already present, to some extent, in the clasaic_work of Infeld and Hull.El That 



was suggested hy some authors several years agcE3E3 and has been shown with some 
detail recently.c3 The ambiguity in the definition of the partner potential is inher- 
ited in the case of Shape Invariance, so one may wonder to what extent it makes 
sense the relation between a potential and its partner characterizing such a kind of 
problems. 

Therefore, two main questions arise. Are there different solutions for the same 
Riccati equation leading to the same partner?. On the other hand, if the Shape 
Invariance condition holds for a certain partner, is it also true for any other possible 
partner?. Our aim is to analyze these questions using, among other thi ngs, t he 
machinery for dealing with Riccati equations developed in previous papers.EiiLJ 

The letter is organized as follows. In the next section we will review briefly the 
concepts of partner potential, when one potential is given, and Shape Invariance. 
In Section 3 we look at the abovementioned ambiguities. In two recent papers Efi 
two alternative factorizations for different choices of the constant d are given. We 
will analyze this point in Section 4, showing that the ambiguity in the factorization 
energy, and hence the existence of some factorizations is due, in turn, to certain 
parameter invariance symmetry of the given potential. Moreover, if the given poten- 
tial can be considered as being part of a pair of Shape Invariant partner potentials, 
more factorizations can be found. We develop these ideas with some illustrative 
examples. 



2. Factorization method and shape invariance 

The simplest way of generating a new exactly solvable Hamiltonian H from a 
known one H is just to consider an invertible bounded operator B, with bounded 
inverse, and defining H = BHB^ 1 . This new Hamiltonian H has the same spec- 
trum as the starting Hamiltonian H. As a generalization^ we will say that two 
Hamiltonian operators H and H are A-related when AH = HA, where A may be 
singular. In this case, if tp is an eigenvector of H corresponding to the eigenvalue E 
and Aip ^ 0, then, at least formally, Aip is also an eigenvector of H corresponding 
to the same eigenvalue E. 
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If wc assume that the intertwining operator A is a first order differential operator, 

A=^- + W(x), and A^ = --^ + W(x), 
dx dx 

then the relation AH = HA, with 



H = -— + V(x), H = -— + V(x), (1) 



£_ ... . ~- d 2 

dx 2 

leads to 

V = -2W' + V, W(V -V) = -W" - V , 

and taking into account the first equation, the second one becomes 2WW = W" + 
V', which can easily be integrated giving 

V = W 2 -W' + d, (2) 

V = W 2 + W' + d, (3) 

where d is an integration constant. The important point here is that H and H, 
given by (Q), can be related by a first order differential operator A of the form given 
above if, and only if, there exists a constant d and a function W such that the pair 
of Riccati equations (||) and (^) are satisfied simultaneously. Moreover, this means 
that both Hamiltonians can be factorized as 

H = A*A + d, H^AA^+d. (4) 

Adding and subtracting equations (^|) and (0) we obtain the equivalent pair 
which relates V and V 

V-d = -{V - d) + 2W 2 , (5) 

V = V + 2W'. (6) 

The function W satisfying these equations is usually called superpotential, the con- 
stant d factorization energy or factorization constant and V and V (resp. H and 
H) are said to be partner potentials (resp. Hamiltonians). 

Gendenshtein took equations (H|) and (|J) as a definition of the functions V, 

V in terms of the function W and some constant d. After, he supposed that W 
did depend on certain set of parameters a, W — W(x,a), and as a consequence 

V = V(x,a) and V = V(x, a) as well. Then, the necessary condition for V(x,a) 
to be of the same form as V(x, a), maybe for a different choice of the values of 
the parameters involved in V, is known as Shape Invariance. More explicitly, it 
amounts to assume the further relation between V(x, a) and V(x, a) 

V(x,a) = V(x,f(a))+R(f(a)), (7) 



where / is an (invertible) transformation on the parameter space a and R is some 
function. The main advantage of these problems is that the complete spectrum of 
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the corresponding Hamiltonians H and H can be found easily.13 Let us remark that 
it is the choice of a and f(a) what defines the different classes of Shape Invariant 
potentials. The function / may be even the identity, f(a) = a.E3 For a more detailed 
information see, for example, Ref. 24 . 

3. On the ambiguity in the definition of the partner potential 

Given one potential function V, the equation (Q) to be solved when searching 
for a superpotential function W, once d is fixed, is a Riccati equation. In general, 
its general solution cannot be found by means of quadratures. However, now we 
only need to compare solutions of the same equation when a particular solution is 
known. In such a case, it is well known that its general solution can be written 
using two quadratures. For a group theoretical explanation of this fact, see, for 
example, Ref. 25. Our aim now is to study a bit further the general solution of (g) 
in that situation and analyze the corresponding possible partner potentials V. 

It is well knownc-K3 that if W p is a particular solution of (^) for some specific 
constant d, the change of variable 

1 1 

with inverse W = W p , (8) 



W p - W v 
transforms (||) into the inhomogeneous first order linear equation for v 

dv 



/lr -2W p v + l, (9) 



which has the general solution 



, s Fexp{2fw p (r,)dT,}dZ + F 
exp{2/ W p (£)d£} 

where F is an integration constant. Therefore, the general solution of (||) reads as 
ii7 i \ wf , e X p{2j x W p (0d^} 

W a (x) = WJx) j . (11) 

fe W {2fw p (r,)dr,}dt + F 

We will review now the concept of partnership given one potential V. We have 
to find a constant d and at least one particular solution W p of the Riccati equation 
(||). Then, "the partner" V is constructed by using (^) or equivalently (Q). But 
these formulas explicitly show that V does depend upon the choice of the particular 
solution of (0) considered. Since the general solution of (||) can be written as 
W g = W p — 1/v, where v is given by (|hJ), the general solution obtained for V g is, 
according to (||), 

d fl 



V ° = V *- 2 d- X [v) ■ (12) 
This answers one of the questions in the introduction: all the partner potentials, 
obtained by using ( |l2| ) are different, apart from the trivial case in which W p and V 
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are constant, because the differential equation (j^) only admits a constant solution 
when W p is constant. 

This implies that "the partner" of one given potential is not a well defined 
concept and it seems better to say that an ordered pair (V, V) is a supersymmetric 
pair of partner potentials if there exists a constant d and a function W such that 
this last is a common solution of the Riccati equations (Q) and (Q) constructed 
with these potentials, respectively. Of course the preceding comment shows that in 
such a case the superpotential function W is essentially unique for each d, which 
moreover makes the problem of ^4-related Hamiltonians be well defined. Note as 
well that this reformulation of partnership comprehends the situation where V is 
the potential we have started this section with, V is one of the functions obtained 
from ( |l2|) for a specific value of the constant F, and W is obtained from (]ll]) for 
the same value of F. 

Now we will show what consequences have this undctcrmination in the subclass 
of Shape Invariant potentials. For that, we should use instead of (J2J) and (||) the 
equations 

V(x,a)-d(a) = W 2 (x, a) - W'(x, a) , (13) 
V(x,a)-d(a) = W 2 (x, a) + W'{x, a) , (14) 

where now the factorization constant depends on the parameter a (see Ref. 24, 
Sec. 3 for details). Consider a particular solution W p (x,a) of equation (|l3|) for 
some specific constant d(a), such that it is also a particular solution of (|l4|), being 
V(x,a) and V(x, a) related by the further condition ([?]). As in the previous case, 
we can consider the general solution of ( jl3| ) starting from W p (x, a), which is 

W g (x,a,F) = W p (x,a)+g(x,a,F), (15) 

where g(x, a, F) is defined by 

( ^ cxp{2f ^ p (C,a)^} 

q{x,a,t) = j , (16) 

rexp{2fw p (r,,a)dr)}d( + F 

being F a integration constant. Note that the particular solution W p (x,a) is ob- 
tained from (15) as F — > oo. Then, inserting W g (x, a, F) into an equation like (fTi)) 



we obtain the general family of partner potentials 

V(x,a,F) = V(x,a)-2g'(x,a,F). (17) 

The question now is whether the condition (Q) is maintained when we consider 
V(x, a, F) and V(x, a) instead of V(x, a) and V(x, a). Then, we ask for 

V(x, a, F) = V(x, f(a)) + R(f(a), F) , (18) 

for some suitable F, where / is the same as in (Q), and R( f(a), F) is a number not 
depending on x, maybe different from the R(f(a)) of (0). Taking into account (fjj) 
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and (p7|), the equation (|l8|) reads as 

2g'(x,a,F) =R(f(a),F) - R(f(a)) , 

that is, 2g'(x, a, F) should be a constant, which we name as k for brevity. Integrating 
respect to x we obtain 2g{x, a, F) = kx + I, being I another constant depending at 
most on a and F. Since g{x, a, F) is, on the other hand, given by (|l6|), it follows 

T < , s i 2 exp{2 f x W„(£,a)<2£) , 

y ex P {2y w 1J (»,,o)d» 7 }de+F= — P1 L + y • ( 19 ) 

Differentiating this last equation, and solving for W p (x, a) we obtain 

1 / 2k \ 

where we have made explicit that the parameter space should be a = {k, I}. Intro- 
ducing this expression in ( |l9[ ) and performing the integrations, we obtain 

— 2 e ~ x (kx+2l)/4 , pi _ _2 e -x(/cx+2;)/4 

and hence F = 0. Now we have to check whether this particular case we have 
found, which is the only candidate for fulfilling (|1S|), satisfy our hypothesis (Q). 
The partner potentials defined by W p (x, k, I) and equations (|l3|), (|lj) are 

fc, I) - d(fc, /) = Wp(x, k, I) - Wp(x, k, I) 

V(x, k, I) - d{k, I) = W^{x, k, I) + W p (x, k, I) 

Now, we have to find out whether there are some transformation of the parameters 
{k, 1} such that the condition ^ be satisfied. Denoting the transformed parameters 
as {ki,l\} for simplicity, we have 

v { x, k, i) - v(x, k lM ) = d( k , i) - d {k M -\-\ (jj^y + (teTo^) 
+^-((k - h)x + 1 - h)((k + h)x + i + h). 

16 

The right hand side of this equation must be a constant and therefore, each of the 
different dependences on x must vanish. The term ((k— k\)x+l— li)((k+ki)x+l+li) 
vanish for the combinations k\ = —k, h = —I or k\ = k, l\ = I, apart form the case 
k\ = —k\ = k = 0, which will be studied separately. However, the term 

3B k 2 
1 

{kix + h) 2 (kx + l) 2 

is equal to Ak 2 /{kx + I) 2 for both combinations and does not vanish. Then, the 
Shape Invariance hypothesis is not satisfied. In the case of k = we have that the 



(kx + Q 2 3fc 2 
16 + A(kx + l) 2 ' 

{kx + l) 2 k 2 k 

16 A{kx + l) 2 ~ 2 ' 
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corresponding W p (x, a) is a constant and hence provides the trivial case where the 
corresponding partner potentials are constant as well. 

This answers the other question of the introduction, and it is closely related with 
the previous one. That is, if the Shape Invariance condition holds for a possible 
partner, then it does not hold for any other choice of partner, apart from the trivial 
case where all the involved functions are constant. 

As a consequence of all the previous, it would be better to reformulate the Shape 
Invariance condition (Q) in terms of appropriate W and d only. Now, considering 



a particular common solution W(x,a) of (13) and ( |14| ) for some d(a), jointly with 
(0) allows to write this last condition as 

W 2 (x, a) - W 2 {x, /(a)) + W'(x, /(a)) + W'(x, a) = R(f(a)) , (20) 

where i?(/(a)) = d(f(a)) — d(a). This way, beginning from W(x, a) and d(a) which 
solve (|2^) for some /, we will obtain through (|l^) and (|l4|) well defined Shape 
Invariant partner potentials (V(x, a), V{x, a)) by construction. In the celebrated 
article by Infeld and Hulffl the key point is indeed to solve an equation of type ( |20| ) 
(see their equation (3.1.2)). Similarly, in Ref. 27, Sec. VI, the main point in the 
classification of Shape Invariant potentials they propose is to find solutions of an 



equation of type ([20|) (see their formula 
a central role in more recent articlesB 



2.22)). An equation of type (|2C| ) also plays 
where some of the results of Infeld and 



Hull are reviewed and put in connection with the concept of Shape Invariance, 
and Shape Invariant partner potentials depending on n parameters transformed 
by translation are found, giving a solution to a previously unsolved problemHj 
However, it seems that a justification of why it is necessary to solve an equation 
of type (pp| ) when searching for well defined Shape Invariant partner potentials has 
not been given explicitly in the literature up to now. 

4. Parameter invariance and Shape Invariance: existence of several fac- 
torizations 

We will analyze in this section what happens if there exists a transformation 
in the parameter space, g : a i— » g(a) such that leaves the potential V(x, a) in 
( |l3|) invariant. Then, whenever (W(x,a), d(a)) is a solution of fll3]), we will have 
another different solution provided W(x, g(a)) ^ W(x,a). In fact, if we transform 
all instances of a in ( |l3| ) by the map g, and use such an invariance property, it 
follows that we have another solution (W(x, g(a)), d(g(a))) of (|l3| ) in addition to 
(W(x,a), d(a)). Inserting each of these pairs into ( fTi]) we will obtain in general 
different partner potentials V(x,g(a)) and V(x,a) of V(x,a). This also gives an 
example of the fact that there may exist several different constants d such that we 
could find a particular solution W of an equation of type (H) or (|l3| ) for a fixed V. 

Another interesting case in which new factorizations can be generated from 
known ones is when we have a pair of partner potentials V(x,a) and V(x,a) sat- 
isfying the Shape Invariance condition (0), properly understood. In this case this 
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condition shows that 

V(x,a) = V(x,f- 1 (a))-R(a), 
or in terms of the Hamiltonians, 

H(a)=H(f- 1 (a))-R(a) 7 

which provides an alternative factorization for H(a): 

H(a) = (-^+W(x,f- 1 (a))^ (~ + W(x, r\a))^ + d^ 1 (a)) - R(a) , 

where it has been used (f|) with A(a) = ^ + W(x, a) and A\a) = — ^ + W(a:, a). 
So, had we started onfo/ with the potential V(af, a) of this paragraph, we would have 
been able to find a factorization of H(a) as a product of type A^(a)A(a) + Const, 
and another as a product A(/ _1 (a))j4J(/ _1 (a)) + Const., being these last constants 
different in general. 

Of course one could have both situations of the preceding paragraphs at the 
same time. We shall illustrate them in the next subsection. 



4.1. Illustrative examples 

As a first example we will explain the four-way factorization of the isotropic har- 
monic oscillator, introduced in Ref. 26, pp. 388-9. In their notation, the potential 
and Hamiltonian of interest are 



1(1 + 1) , _ 2 ufn _ d 2 

dr 2 



y(r,0 = -Hr^ + ^> H(l) = -— 5 + V(r,l). 



where the independent variable is r E (0, oo) and the set of parameters is simply I. 
Their factorization (6) is 

H(l) = (~+ l -+r)(±+ l - + A-(2l-l), (21) 



dr r J \dr r 

from where it is suggested that W(r, I) = ~ + r. Substituting it in V(r, I) = 
W 2 (r,l) - + d(l) we obtain d(l) = -(21 - 1), so © is the appropriate 

version of our (Q) as expected. Now, as the potential V(r, I) is invariant under 
the map g : I i— v — I — 1, we will obtain a new solution (W(r,g(l)) ,d(g(l))) — 
(W(r, — I — 1) , d(— I — 1)) of the equation 

V(r,l) = W 2 ~^+d. 

dr 

But W(r, g(l)) = W(r, -l-l) = -^±+r and d(g(l)) = d(-l - 1) = 21 + 3, which 
is exactly what corresponds to the factorization (4) of Ref. 26. The factorizations 
(5) and (7) loc. cit. are related in a similar way; (7) is obtained from (5) by means 
of the change g : I ^ — I — 1 as well. 
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As far as the relation between their factorizations (6) and (5) loc. cit. is con- 
cerned, we have already seen that, from their factorization (6), here reproduced as 
(^H), it follows W(r, I) = £ +r, and thus, the corresponding V(r, I) through ( |l4| ) is 

V(r, I) = W*(r, I) + + d{l) = l l_R + r2 + 2 . 

ar r z 

Then it is very easy to check that V(r, I) = V(r, /(/)) + R(f(l)), where R(l) = 2 for 
all I, and / is defined either by f(l) = I — 1 or / (Z) = —I. We obtain 

iT(Z) = H(l + 1) - fl(Z) , V(r, /) - t/(r, I + 1) - , 

and 

H(l) - #(-/) - R(l) , K(r, = V(r, -I) - R(l) , 

as well. In this way the factorization (5) of Ref. 26 is achieved. 

As a second example we will consider the modified Pdschl-Teller potential, an- 
alyzed in an interesting recent article.Elj The potential is now 

V(x, a) X) = -a 2 X{X ~ 1} , (22) 
cosh ax 

where x £ (— oo, oo) and a > 0, A > 1 are two real parameters. 

Two different particular solutions (W(x, a, A), d(a, A)) of the Riccati equation 

W 2 - W = V(x, a, A) -d, 

have been found in Ref. 20, p. 8450, namely, 

(Wi(x, a, A), di(a, A)) = (— A a tanh 2 ax, — \ 2 a 2 ) , 
(W 2 (x,a,X),d 2 (a,X)) = (-(1 - A) a tanh 2 ax, -(1 - A) 2 a 2 ) . 

It is clear that the second pair is obtained from the first by means of the parameter 
transformation g : (a, A) i— > (a, 1 — A). The reason is that V{x,a,X) is invariant 
under g, or more precisely, its factor A(l — A). 

The associated partner potentials V(x, a, A) obtained using (|l4|), are 

VifoajA) = W?(g,a,A) +W[(g 1 a,A) +di(a,A) = -a 2 A ^ ^ 



cosh^ 



y 2 (a;,a,A) = W 2 {x, a, A) + Wa(a;, a, A) + e? 2 (a, A) = -a 2 — — ^ — ^. 

cosh aa; 

We see that both of the previous functions are just second degree monic polynomials 
in A, with roots spaced one unit, times —a 2 / cosh 2 ax, like V(x, a, A) itself. It is 
then obvious that a translation of the type Ah- >A — 6orAi— > c — A should transform 
Vi(x, a, A) and Vi{x,a,\) into V{x,a,X). This is in fact so, since V{x,a,X) — 
V\(x, a, / _1 (A)), where / is defined either by /(A) = A — 1 or /(A) = —A, and 
similarly V(x, a, A) = V 2 (x, a, / _1 (A)) when /(A) = A - 1 or /(A) = 2 - A. 
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In this way one could propose other different factorizations for the potential 
V(x,a, A), being able in principle to make a differential operator analysis for this 
potential similar to what it is done in Ref. 26 for the first example of this Subsection. 
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